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ABSTRACT 

This  tutorial  paper  on  polarimetnc  definitions  uses  N-port  network  theory  to  prove  results  that 
are  frequently  stated  without  proof  in  the  literature,  including  the  symmetry  of  the  polarization 
scattering  matrix  in  backscattering,  the  equation  for  received  voltage,  and  the  equations  defining 
Stokes  vectors  and  backscatter  Mueller  matrices.  Appropriate  N-port  networks  are  defined 
for  a  single-polarized  antenna,  a  dual-polarized  antenna,  and  a  backscattering  target.  An 
important  result  is  demonstrated:  reciprocity  (symmetry  of  the  polarization  scattering  matrix) 
is  meaningful  only  in  the  context  of  the  monostatic  convention,  for  which  the  coordinate  system 
is  the  same  for  both  transmit  and  receive.  This,  in  turn,  implies  a  change  in  the  handedness 
of  the  coordinate  system,  so  that  scattered  fields  and  Stokes  vectors  are  expressed  in  opposite 
handedness  from  incident  fields  and  Stokes  vectors. 

These  results  are  used  as  tools  in  the  derivation  of  a  new  technique  for  measuring  depolar¬ 
ization  in  foliage  penetration  (FOPEN).  Four  polarimetric  active  radar  calibrators  (PARCs)  are 
used  with  an  algorithm  to  completely  measure  the  Mueller  matrices  describing  downward-going 
and  upward-going  foliage  penetration. 
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1.  INTRODUCTION  AND  CONCLUSIONS 


This  is  a  tutorial  paper  on  polarimetric  definitions  (scattering  matrices,  Stokes  vectors,  and 
backscatter  Mueller  matrices).  The  paper  contains  a  number  of  results  which  have  proved 
useful  to  the  writer,  and  also  introduces  the  notion  of  coordinate-system  handedness,  which,  it 
is  felt,  is  not  adequately  treated  in  the  literature. 

These  results  are  used  to  develop  a  new  technique  for  measuring  depolarization  in  fo¬ 
liage  penetration  (FOPEN):  four  polarimetric  active  radar  calibrators  (PARCs)  are  used  with 
an  algorithm  to  completely  measure  the  Mueller  matrices  describing  depolarization  during 
downward-going  and  upward-going  foliage  penetration. 

The  major  results  derived  in  this  paper  include  the  following: 

1 .  The  polarization  scattering  matrix  for  backscatter  is  symmetric  only  if  the 
incident  field  and  backscattered  field  are  expressed  in  coordinate  systems 
that  have  opposite  handedness.  For  example,  if  the  transmit  field  is  ex¬ 
pressed  in  a  right-handed  system,  the  backscattered  field  must  be  expressed 
in  a  left-handed  system.  If  one  transforms  the  backscattered  field  so  that  it 
also  is  in  a  right-handed  system,  then  the  scattering  matrix  is  not  symmetric. 

2.  The  voltage  measured  by  a  receive  antenna  can  be  written  as  a  simple 
dot  product  rTe  between  the  receive  antenna  complex  height  r  and  the 
backscattered  field  e,  where  e  is  expressed  in  a  left-handed  coordinate 
system.  This  is  not  an  inner  product  (e.g.  r}e),  where  the  superscript 
dagger  denotes  conjugate  transpose,  but  rather  a  simple  dot  product. 

3.  The  receive  distortion  matrix  B  of  a  dual-polarized  antenna  combined  with 
an  orthomode  transducer  is  equal  to  the  transmit  distortion  matrix  C.  This 
result  follows  from  reciprocity,  and  is  therefore  not  generally  true  for  a 
complete  receiver,  which  may  include  nonreciprocal  elements  such  as  cir¬ 
culators. 

4.  The  equation  for  the  receive  antenna's  measured  power  can  be  written  in 
terms  of  the  Stokes  vectors  of  the  receive  antenna  and  the  backscattered 
field.  The  backscattered-field  Stokes  vector  must  be  expressed  in  a  left- 
handed  coordinate  system. 

5.  The  backscatter  Mueller  matrix  (also  known  as  the  Stokes  reflection  matrix) 
is  symmetric  only  if  the  incident  and  scattered  Stokes  vectors  are  expressed 
in  coordinate  systems  with  opposite  handedness. 
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6.  The  transmission  Mueller  matrix  (for  example,  the  matrix  describing  depo¬ 
larization  through  foliage)  is  not  symmetric,  and  therefore  has  sixteen  and 
not  nine  independent  elements 

Section  2  uses  N-port  network  theory  to  derive  several  simple  relationships  between  complex 
two-vectors  and  polarization  scattering  matrices.  Section  3  applies  the  results  from  Section  2 
to  demonstrate  relationships  between  Stokes  vectors  and  backscatter  Mueller  matrices.  Section 
4  extends  the  backscatter  results  to  forward  scatter,  and  then  uses  the  results  as  tools  to  develop 
a  technique  for  measuring  depolarization  in  foliage  penetration. 
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2.  N-PORT  NETWORK  THEORY  APPLIED  TO 
BACKSCATTERING  AND  ANTENNAS 


Engineering  progress  can  frequently  be  made  by  exploiting  analogies  with  results  developed 
in  similar  fields.  Such  an  analogy  is  proposed  here,  between  radar  antennas  and  scattering 
objects  on  the  one  hand,  and  a  branch  of  linear  systems  theory  called  N-port  network  theory 
on  the  other  hand.  In  this  section  we  will  model  a  backscattering  object  as  a  two-port  network, 
a  single-polarized  antenna  as  a  three-port  network,  and  a  dual -polarized  antenna  as  a  four-port 
network.  In  order  to  do  this,  some  basic  results  from  N-port  network  theory  ([2],  e.g.)  are  first 
recalled. 


An  N-port  network  (see  Figure  2-1)  can  be  characterized  by  a  scattering  matrix,  which 
relates  backscattered  to  incident  waves.  At  port  m,  the  total  electric  field  can  be  expressed 
as  Em  =  (am  +  6m)em,  while  the  total  magnetic  field  is  Hm  =  (a™  -  bm)hm,  where  em 
and  hm  are  the  unit  vectors  for  the  electric  and  magnetic  fields  at  the  m-th  port,  a™  is  the 
complex  coefficient  of  the  inward  wave  at  the  m-th  port,  and  bm  is  the  complex  coefficient 
of  the  outward  wave  at  the  m-th  port.  Note  that  both  inward  and  outward  waves  are  defined 
with  respect  to  the  same  coordinate  system  (em,hm).  The  scattering  matrix  5  relates  these 
coefficients  as  follows: 


h  =  Sfl 


(2.1) 
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a  i 

where  b  = 

■ 

and  a  - 

• 

aN 

and  where  the  element  is  the  amplitude  of  an  outward  wave  at  port  i  resulting  from  a  unit 
inward  wave  at  port  j,  with  all  other  ports  “matched”  (i.e.  with  no  other  inward  waves).  In 
this  expression  a  and  b  are,  in  general,  complex  N-vectors. 

An  important  result  from  N-port  network  theory  is  that  when  an  N-port  network  is  passive 
and  reciprocal,  its  scattering  matrix  5  is  symmetric,  i.e.,  S,:  =  S^.  This  result  will  be  used 
extensively  in  this  paper, 

2.1  The  Backscatterer  as  a  Two-port  Network 

N-port  network  theory  can  be  used  to  model  a  backscattering  object  as  a  two-port  network. 
Define  Ports  1  and  2  as  corresponding  to  horizontally  and  vertically  polarized  electric  fields, 
respectively;  then  the  unit  vectors  ex  and  e2  correspond  to  xt  and  yf,  respectively,  in  Figure  2-2. 
The  particular  choice  for  ex  and  e2  is  not  important;  what  is  important  is  that  the  same 
coordinate  system  is  used  for  both  inward  and  outward  waves  (i.e.  for  the  incident  and 
backscattered  fields).  Since  the  propagation  vectors  for  the  incident  and  backscattered  fields  are 
in  opposite  directions,  this  means  that  the  incident  and  scattered  coordinate  systems  must  have 
opposite  handedness.  We  arbitrarily  choose  a  right-handed  triplet  (£t,  yt,zt)  for  the  transmit 
coordinate  system  and  a  left-handed  triplet  (£,,  y,,  £,)  for  the  backscattered  coordinate  system, 
as  shown  in  Figure  2-2. 

With  these  definitions,  the  polarization  scattering  matrix  A  for  a  backscattering  object  can 
be  viewed  as  the  N-port  scattering  matrix  S  for  the  case  N  =  2.  We  will  now  find  analogies 
for  each  of  the  parts  (q ,,  5,  and  £)  of  Equation  2.1: 

1.  The  vector  a  of  incoming  waves  corresponds  to  the  incident  field  two- 
vector  t,  which  we  use  to  denote  the  field  transmitted  by  the  transmit 
antenna.  Note  that  i  is  expressed  in  a  right-handed  coordinate  system. 

2.  The  vector  b  of  the  outgoing  waves  corresponds  to  the  backscattered  wave 
£,  (the  tilde  denotes  that  £  is  expressed  in  a  left-handed  coordinate  system). 

3.  Finally,  the  N-port  scattering  matrix  S  corresponds  to  the  polarization  scat¬ 
tering  matrix  A  for  the  backscattering  target. 

Combining  these  substitutions  with  the  original  N-port  equation  (Equation  2.1),  we  get: 
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ys  =  y  i 


2S=  Z. 


SCATTERED 

COORDINATE 

SYSTEM 

[LEFT-HANDED] 


TRANSMIT  (INCIDENT) 
COORDINATE  SYSTEM 
[RIGHT-HANDED] 


Figure  2-2.  Geometry  for  a  seatterer  as  a  two-port  network. 

£  =  At  (2.3) 

As  an  example  of  the  above  result,  consider  a  left -circular  wave  (expressed  in  a  right-handed 
coordinate  system)  which  is  assumed  to  be  incident  on  a  sphere.  The  incident  wave  can  be 

expressed  as  t  —  4-  *  ,  and  the  scattering  matrix  of  a  sphere  is  S  =  1  ?  ;  thus  the 

j  0  1 

scattered  field  is  e  =  At  =  \  .  Physically,  we  know  that  the  scattered  field  must  be 


right-circular  polarized.  Thus,  the  same  complex  two-vector  4-  refers  to  left-circular  in 

J 

the  right-hai  ded  coordinate  system  and  right-circular  in  the  left-handed  coordinate  system. 


2.2  The  Single-polarized  Receive  Antenna  as  a  Three-port  Network 

N-port  network  theory  can  also  be  applied  to  modelling  an  antenna.  In  this  section,  we  will 
model  a  single-polarized  antenna  as  a  three-port  network.  Define  Ports  1  and  2  as  corresponding 
to  the  horizontally  and  vertically  polarized  electric  fields,  respectively,  transmitted  or  received 
by  the  antenna.  Also,  define  Port  3  as  the  waveguide  input/output  at  which  the  antenna’s 
complex  voltage  is  defined. 

Assume  that  the  antenna  in  question  is  ofil  fitted  out  with  an  orthomode  transducer.  (If  it 
were,  it  would  have  to  be  modelled  as  a  four-port  network;  this  case  will  be  treated  in  Section 
2.3.) 
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Let  us  assume  that  the  the  unit  vectors  ex  and  e2  correspond  to  xt  and  y,,  respectively, 
in  Figure  2-3.  This  means  that  transmitted  waves  are  expressed  in  a  right-handed  coordinate 
system  and  received  waves  are  expressed  in  a  left-handed  coordinate  system. 

PORT  2 


y : 


Figure  2-3.  Geometry  for  a  single-polarized  antenna  as  a  three-port  network. 


Consider  the  antenna  operating  as  a  transmitter.  Assume  there  is  a  unit  incoming  wave 
at  Port  3  (i.e.,  the  signal  to  be  transmitted  coming  into  the  antenna  waveguide  input/output). 
Assume  also  that  there  are  no  waves  incident  on  the  antenna  from  the  rest  of  the  world,  so 
that  the  incoming  waves  at  Ports  1  and  2  are  zero.  The  g,- vector  (the  complex  three-vector 
describing  the  coefficients  of  the  incoming  waves  to  the  three  ports  of  the  antenna)  can  therefore 
be  expressed  as  [  0  0  1  ]T,  so  that  (according  to  Equation  2.1)  the  resulting  6-vector  (the 
complex  three-vector  of  outgoing  coefficients)  is  the  last  column  of  the  5  matrix: 


6  = 


(2.4) 


But  the  first  two  elements  of  h  are  the  horizontal  and  vertical  components  of  the  field  that 
this  antenna  transmits.  Call  these  gx  and  y2.  (The  field  £  is  a  dummy  parameter;  it  is  set  to 
i  if  the  antenna  is  to  be  used  as  a  transmit  antenna,  and  to  r  if  it  is  to  be  used  as  a  receive 
antenna.)  Thus,  Sx3  —  gx  and  S23  =  y2.  Because  the  5  matrix  is  symmetric,  we  also  have 
=  gx  and  S32  -  g2. 
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Now  consider  the  operation  of  the  antenna  as  a  receiver.  Here  the  analogy  with  the  N-port 
network  is  as  follows:  the  role  of  the  incoming  wave  is  played  by  the  electric  field  incident 
on  the  antenna  from  the  outside  world,  while  the  role  of  the  outgoing  wave  is  played  by  the 
signal  sent  down  the  waveguide  from  the  antenna  to  the  receiver. 

Assume  that  an  incident  wave  with  horizontal  component  e1  and  vertical  component  e2 
arrives  at  Ports  1  and  2,  and  that  no  signals  are  coming  into  Port  3  (the  antenna  waveguide 
input/output)  from  the  receiver.  Thus  the  incoming  wave  vector  to  be  used  with  the  three-port 
network  is  a  =  [  e,  e2  0  ]T.  Let  us  examine  the  value  of  the  outgoing  wave  in  Port  3,  i.e., 
the  antenna  received  voltage.  This  is  the  third  component  of  the  b- vector,  which  is  equal  to 


V  —  b3  —  S3icii  -1-  532a2  -h  53303  --  ffjej  4-  g2e2 


(2.5) 


This  result  can  then  be  applied  to  the  receive  antenna  (for  which  g  is  taken  to  be  r): 


V  —  rTe 

Note  that  this  is  not  an  inner  product;  i.e.  neither  r  nor  e  is  to  be  conjugated. 


(2.6) 


2.3  The  Dual-polarized  Antenna  as  a  Four-port  Network 

An  orthomode  transducer  is  a  microwave  device  with  an  input/output  that  can  support  elec¬ 
tric  fields  in  either  or  both  of  two  orthogonal  modes  (these  modes  are  typically  referred  to  as 
horizontally  and  vertically  polarized),  as  well  as  two  waveguide  connections.  A  dual-polarized 
antenna  is  just  an  orthomode  transducer  mounted  onto  a  feed/aperture  structure  capable  of  sup¬ 
porting  both  horizontal  and  vertical  polarizations.  Since  we  count  one  port  for  each  waveguide 
connection  and  one  port  each  for  the  horizontal  and  vertical  components  of  the  radiated  electric 
field,  the  dual-polarized  antenna  can  be  modelled  as  a  four-port  network. 

The  analysis  of  this  four-port  network  will  closely  parallel  that  in  Section  2.2.  Define  Ports  1 
and  2,  respectively,  as  corresponding  to  the  horizontally  and  vertically  polarized  electric  fields 
transmitted  or  received  by  the  antenna.  Define  Ports  3  and  4,  respectively,  as  the  horizontally 
and  vertically  polarized  output  of  the  orthomode  transducer. 

Let  us  assume  that  the  the  unit  vectors  and  e2  correspond  to  xt  and  yt,  respectively, 
in  Figure  2-4.  This  means  that  transmitted  waves  are  expressed  in  a  right-handed  coordinate 
system  and  received  waves  are  expressed  in  a  ieft-handed  coordinate  system. 

Consider  the  antenna  operating  as  a  transmitter.  In  this  mode,  the  two  waveguide  connections 
(horizontal  and  vertical)  of  the  orthomode  transducer  are  used  as  inputs.  Although  signals  can 
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PORT  2 


y  i 


Figure  2-4.  Geometry  for  a  dual-polarized  antenna  as  a  four-port  network. 

be  present  in  both  waveguides  at  the  same  time,  we  will  analyze  the  cases  for  which  one  one 
or  the  other  waveguide  input  is  used,  and  then  appeal  to  linear  superposition  to  find  the  field 
radiated  when  both  waveguide  inputs  are  used  simultaneously. 

First  we  assume  that  only  the  horizontal  waveguide  input  is  used. 

Assume  there  is  a  unit  incoming  wave  in  Port  3  (i.e.,  assume  that  the  signal  to  be  transmitted 
is  coming  into  the  orthomode  transducer  horizontal  input).  Assume  also  that  there  are  no  waves 
incident  on  the  antenna  from  the  rest  of  the  world,  so  that  the  incoming  waves  in  Ports  1,2,  and 
4  are  zero.  The  a-vector  (the  complex  four-vector  describing  the  coefficients  of  the  incoming 
waves  to  the  four  ports  of  the  antenna)  is  thus  [0  0  1  0  ]r,  so  that  the  resulting  6-vector 
(the  complex  four-vector  of  outgoing  coefficients)  is  (by  Equation  2.1)  equal  to  the  third  column 
of  the  5  matrix: 


h 


Sq.  = 


S 13 

S23 

£33 

S43 


(2.7) 


But  the  first  two  elements  of  h  are  the  Port  1  and  Port  2  outputs,  i.e.,  the  horizontal  and 
vertical  components  of  the  field  the  antenna  transmits  when  the  nominal  transmit  polarization 
is  horizontal.  Call  these  two  elements  Cn  and  C 21,  using  the  notation  of  [4].  The  electric  field 
actually  transmitted  when  we  input  only  to  the  horizontal  waveguide  is  £1  =  (  C\  1  C21  ]r- 
Thus  S\3  —  C\i  and  Sj 3  =  C21  • 
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Next,  we  assume  that  only  the  vertical  waveguide  input  is  used. 

The  analysis  is  similar  to  that  used  for  the  horizontal  waveguide.  In  this  case,  only  the  vertical 
input  of  the  orthomode  transducer  is  stimulated,  so  the  a-vector  is  equal  to  [  0  0  0  1  ]r. 
The  resulting  6- vector  is  (by  Equation  2.1)  equal  to  the  fourth  column  of  the  S  matrix: 


b  = 


Sa  = 


S14 

S24 

S34 

S44 


(2.8) 


Thus,  using  the  notation  of  [4],  when  we  input  to  the  vertical  waveguide  (i.e.  into  Port  4), 
the  actual  transmitted  field  is  C2  =  [  Cu  C22  ]r;  since  these  two  coefficients  are  equivalent 
to  the  first  two  elements  of  the  fourth  column  of  S ,  we  have  =  Ci2  and  524  =  C22. 

What  has  been  shown  so  far  is  that  the  upper-right  two-by-two  submatrix  of  the  four-by-four 
scattering  matrix  is  the  transmit  distortion  matrix  C  =  [£Zi£22]  defined  in  [4], 

Now  consider  the  operation  of  the  antenna  as  a  receiver.  Assume  that  an  incident  wave 
with  horizontal  component  el  and  vertical  component  e2  arrives  at  Ports  1  and  2  respectively, 
and  that  no  signals  are  coming  into  Ports  3  and  4  (the  orthomode  transducer  waveguide  ports) 
from  the  receiver.  Thus  the  incoming  wave  vector  to  be  used  with  the  four-port  network 
analogy  is  a  =  (  ex  e2  0  0  ]T.  Let  us  examine  the  value  of  the  outgoing  wave  in  Port  3, 
i.e.,  the  antenna  output  voltage  in  the  horizontal  orthomode  transducer  port.  This  is  the  third 
component  of  the  b- vector,  which,  using  the  symmetry  of  the  5  matrix,  is  given  by: 


V  =  63 


<^31al  +  S32CL2  +  S33<l3  +  ‘§’3404 

5i3ai  4-  S23CL2  +  S330.3  +  £4304 


—  Cllel  +  C21e2 

=  aU  (2-9) 

Similarly,  the  output  voltage  in  Port  4  (vertical)  is  V  =  64  =  Q^e. 

The  complex  numbers  63  and  64  are  the  horizontal  and  vertical  components  of  the  distorted 
receive  field  (where  the  term  “distorted”  is  the  tern  used  in  [4]  to  refer  to  a  measurement 
which  differs  from  its  nominal  value).  If  63  and  64  are  arranged  as  a  complex  two-vector, 
which  we  call  £',  then 


-/ 


£  = 


(2.10) 
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But  the  receive  distortion  matrix  B  is  defined  in  [4]  by  the  equation  e  = 
CTe  =  BTe  for  arbitrary  choice  of  e,  which  in  turn  implies  that  B  =  C. 


63 

64 


BTe.  Thus 


We  have  now  shown  that  a  dual-polarized  antenna  has  B  =  C\  i.e.,  the  transmit  and  receive 
distortion  matrices  are  equal.  This  result  is  not  generally  true  if  the  analysis  is  extended  beyond 
the  antenna  to  include  nonreciprocal  elements  such  as  circulators  in  the  receiver  chain. 
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3.  THE  FLIP  MATRIX,  STOKES  VECTORS,  AND  THE 

MUELLER  MATRIX 


In  this  section  the  characterization  of  scatterers  and  antennas  in  terms  of  complex  two-vectors 
developed  in  Section  2  will  be  used  to  derive  simple,  explicit  definitions  for  Stokes  vectors  and 
backscatter  Mueller  matrices  in  terms  of  Kronecker  products.  Given  these  simple  expressions, 
we  will  quickly  derive  a  number  of  results  frequently  stated  without  proof: 

1 .  The  received  voltage  expressed  in  terms  of  Stokes  vectors  for  the  receive 
antenna  and  backscattered  field. 

2.  The  backscattered  Stokes  vector  expressed  in  terms  of  the  backscatter 
Mueller  matrix. 

3.  The  symmetry  of  the  backscatter  Mueller  matrix. 

4.  The  “trace  rule”  for  backscatter  Mueller  matrices. 


3.1  The  Flip  Matrix 

We  begin  with  the  monostatic  convention  shown  in  Figure  2-2,  with  a  right-handed  triplet 
(it,  yt,  zt)  for  the  transmit  coordinate  system  and  a  left-handed  triplet  (£,,  y,,  i.)  for  the 
scatter  coordinate  system.  It  was  shown  in  Section  2  that  the  backscattered  field  £  expressed 
in  the  left-handed  coordinate  system  and  the  received  voltage  V  are  given  by: 


£  —  At 

V  =  rT  £  (3.1) 

where  t  and  r  are  expressed  in  the  right-handed  transmit  coordinate  system,  and  £  is  expressed 
in  the  left-handed  backscatter  coordinate  system  (see  Figure  3-1).  (As  in  Section  2,  the  tilde 
is  used  to  show  that  a  field  is  expressed  in  a  left-handed  system.) 

Workers  in  the  field  usually  use  scattered  Stokes  vectors  expressed  in  a  right-handed  coordi¬ 
nate  system.  If  one  wishes  to  change  the  basis  of  the  scattered  field  to  a  right-handed  system, 
(-£,,  y,,  z,)  one  can  flip  the  x-axis  by  premultiplying  by  a  flip  matrix  F: 


£  =  Fi 


(3.2) 


where  F 


-1  0 

0  1 


(3.3) 


11 
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Note  that  while  the  matrix  A  is  symmetric,  the  matrix  FA  is  not. 

The  next  section  demonstrates  how  to  write  Stokes  vectors  for  both  transmitted  and  backscat- 
tered  fields. 


3.2  Stokes  Vectors 

The  transmit  and  receive  Stokes  vectors  (i.e.  the  Stokes  vectors  associated  with  the  fields 
that  would  be  transmitted  by  the  transmit  and  receive  antennas)  are  [1]: 


where 


st  =  T(t®t") 

Sr  =  T(t  ®  r") 

10  0  1 
T=  1  0  0-1 
0  110 
0  j  -j  0 


(3.4) 

(3.5) 


(3.6) 


and  where  the  symbol  <g>  denotes  the  direct  or  Kronecker  product  [5,  6].  We  use  the  left 
Kronecker  product:  if  A  and  B  are  two-by-two  matrices,  then  the  Kronecker  product  of  A  and 
B  is 


A  ®  B 


CLliB  Q12B 
CL2\B  O-22B 


(3.7) 


Note  that  the  matrix  T  is  unitary,  since  T_1  = 

The  scattered  Stokes  vector  (expressed  in  a  left-handed  coordinate  system)  is  given  by  (see 
the  Appendix  and  [3]) 


ie  = 

=  T‘(FFe)  ®  (F*F‘e‘) 

=  r(F®r)(Fi«rf) 

=  r*(F®F*)7’-1r(£®£*)  (3.8) 
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(Here  we  have  used  the  fact  that  the  matrix  product  FF  is  equal  to  the  identity  matrix.) 

Equation  3.8  can  be  expressed  more  compactly  by  ^  =  Z^,  where  ^  =  T(e  ®  e")  is 
the  backscattered  Stokes  vector  expressed  in  a  right-handed  coordinate  system,  and  where  the 
4-by-4  real  matrix  Z  is,  as  will  be  shown  in  Equation  3.12,  the  Mueller-matrix  equivalent  of 
the  flip  operation: 


Z  =  T'(F  ®  F*)r_1 


10  0  0 

0  10  0 

0  0-1  0 

0  0  0  -1 


(3.9) 


3.3  Received  Power  and  The  Backscatter  Mueller  Matrix 

Equations  3.4  and  3.8  can  be  combined  to  show  the  relationship  between  received  power, 
the  backscattered  Stokes  vector,  and  the  backscatter  Mueller  matrix.  The  received  power  is 
given  by 


P=  \V\2  = 


(rT|)  ®  (rfe') 

(lt  ®  £’) 

(r  ®  L*)T(T‘)~lT'(e  ®  e‘) 

i[(r  ®r*)rTT]ie 
i(Tr  ®  r*)TIe 


and  the  scattered  Stokes  vector  is  given  by 


ie  =  T'i  ®  £* 

=  T-(Ai®AT) 

=  T'{A  ®  A*)(t  ®r) 

=  (TM®  A'T-l)(TL®f) 

=  drA®A-T')i, 


(3.10) 


(3.11) 
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This  latter  result  can  be  written  a si,  =  Mst,  where 


M  =  T’(A®  A')T~l  =  A')T'  (3.12) 

is  the  Mueller  matrix  relating  backscattered  to  incident  Stokes  vectors.  Note  that  the  scattered 
Stokes  vector  is  expressed  in  a  left-handed  coordinate  system,  whereas  the  incident  Stokes 
vector  is  expressed  in  a  right-handed  coordinate  system. 

Combining  Equations  3.10  and  3.12: 


P  =  \i TMl,  (3.13) 

where  both  the  transmit  and  receive  Stokes  vectors  (i.e.,  the  Stokes  vectors  corresponding  to 
both  antennas,  transmit  and  receive)  are  expressed  in  a  right-handed  coordinate  system. 

We  have  developed  several  tools  useful  for  working  with  Stokes  vectors  and  Mueller  matri¬ 
ces,  and  have  already  used  these  tools  to  derive  results  for  the  scattered  Stokes  vector  (Equation 
3.8)  and  the  received  power  in  terms  of  the  Mueller  matrix  (Equation  3.12).  With  these  tools, 
it  is  easy  to  derive  other  results  that  are  frequently  stated  without  proof.  For  example,  let  us 
derive  the  symmetry  of  the  backscatter  Mueller  matrix  and  prove  the  “trace  rule”. 

Symmetry  of  the  Backscatter  Mueller  Matrix 

We  can  show  that  the  backscatter  Mueller  matrix  M  is  symmetric,  using  the  symmetry  of 
A  and  the  identities  (M  ®  N)T  =  MT  ®  NT  and  ( MN)T  =  NT MT: 


Mt  =  ^T‘{A®  A')tT] 

=  \t'{At  ®  A')T* 

=  \t’{A®A‘)T' 

=  M  (3.14) 

Since  the  backscatter  Mueller  matrix  is  real  and  symmetric,  it  can  have  at  most  ten  indepen¬ 
dent  parameters,  since  the  six  elements  below  the  main  diagonal  are  equal  to  their  counterparts 
above  the  main  diagonal.  We  next  show  that  the  number  of  independent  elements  is  nine  or 
fewer,  because  the  elements  on  the  main  diagonal  are  not  independent  of  each  other. 


Let  us  prove  that  the  (1,1)  element  of  the  backscatter  Mueller  backscatter  matrix  M  is  equal 
to  the  sum  of  the  other  three  diagonal  elements.  This  relationship  is  called  the  “trace  rule” 
(the  trace  operator  is  the  sum  of  the  diagonal  elements  of  a  matrix)  because  the  relationship 
can  be  stated  using  the  trace  as  follows:  Tr(RM)  =  0,  where  the  matrix  R  is  defined  as 


R 


10  0  0 
0-100 
0  0-10 
0  0  0  -1 


The  proof  involves  simple  substitution  and  calculation.  First  the  substitution: 


Tr(RM)  =  Tr(RT‘A  <g>  AmT~x) 

=  Tr(T~x  RT’  A  ®  A")  (3.15) 

where  the  last  step  follows  from  the  fact  that  Tr(AB)  =  Tr(BA).  Next,  the  value  of  the 
matrix  T~1RT‘  can  be  found  by  direct  calculation: 


T~'RT’  = 


0  0  0  1 
0  0-10 
0-1  0  0 

10  0  0 


This  matrix  can  be  represented  as  a  direct  (Kronecker)  product  of  a  simple  2-by-2  matrix  G 
with  its  own  conjugate: 


T-'RT' 
where  G 

Thus, 

Tr{RM) 


G®G’ 

'  0  -1  ' 

1  0 

7Y[((7®G*)(i4®  A')] 

Tr[{GA)®{G‘A-)} 

Tr(GA)Tr(G'A‘) 
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(3.16) 


=  iTr(GA)f 
—  1^21  ~  -^4.12 12 
=  0 

where  the  first  two  steps  follow  from  standard  properties  of  direct  products  [5]:  (A  ®  B)(C  ® 
D)  —  AC<S>BD,  and  Tr(A®B )  =  Tr(A)Tr(B),  and  the  last  step  follows  from  the  symmetry 
of  A. 

This  completes  the  proof  that  the  trace  is  zero. 
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4.  MEASUREMENT  OF  FORWARD-SCATTER 
DEPOLARIZATION 


The  previous  sections  of  this  paper  described  backscatter  polarimetry,  for  which  the  Mueller 
matrix  reduces  to  a  special  form  known  as  the  Stokes  reflection  matrix.  This  Stokes  reflection 
matrix  is  symmetric  and  obeys  the  “trace  rule”  (the  (1,1)  element  is  equal  to  the  sum  of 
the  other  three  diagonal  elements).  As  a  consequence,  the  matrix  has  not  sixteen  but  nine 
independent  elements. 

In  contrast,  the  Mueller  matrix  that  describes  forward  scatter,  which  we  refer  to  here  as  the 
transmission  Mueller  matrix,  has  all  sixteen  elements  independent. 

Analysis  of  forward  scatter  is  of  paramount  importance  in  foliage  penetration  (FOPEN) 
systems,  in  which  a  radar  wave  is  to  pass  through  tree  foliage,  bounce  off  a  target  hiding  under 
the  trees,  and  pass  back  up  through  the  foliage.  Measurements  of  foliage  depolarization  that 
have  been  described  in  the  literature  do  not  include  a  systematic  mathematical  structure  with 
which  to  evaluate  the  measurements.  Using  the  N-port  and  Stokes-vector  theory  developed 
in  previous  sections  of  this  paper,  we  will  not  only  describe  forward-scatter  depolarization, 
but  also  develop  an  algorithm  for  systematically  measuring  depolarization  using  polarimetric 
active  radar  calibrators  (PARCs). 

Section  4.1  describes  the  problem  which  is  to  be  solved.  Section  4.2  develops  the  tools 
needed  for  describing  forward-scatter  depolarization:  Section  4.2.1  derives  the  Stokes  reflec¬ 
tion  matrix  for  a  PARC,  and  Section  4.2.2  demonstrates  the  equality  of  the  Mueller  matrices 
describing  upward-going  and  downward-going  foliage  penetration.  Section  4.3  presents  a  new 
depolarization  measurement  algorithm,  with  the  relative  (depolarization)  part  in  Section  4.3.1 
and  the  absolute  (attenuation)  part  in  Section  4.3.2. 

4.1  Problem  Definition 

Assume  that  the  polarization-changing  properties  of  foliage  can  be  summarized  by  a  trans¬ 
mission  Mueller  matrix  N  which  describes  how  a  Stokes  vector  ^  incident  on  the  foliage  is 
changed  to  the  transmitted  Stokes  vector  s.  which  emerges  on  the  other  side  of  the  foliage: 


s2  =  Nsx  (4.1) 

Imagine  a  situation  in  which  a  radar  wave  originates  above  a  treeline  in  an  airborne  trans¬ 
mitter,  propagates  through  the  trees,  bounces  off  a  point  target  under  the  trees,  and  propagates 
back  up  through  the  trees  and  then  back  to  the  radar.  Denote  the  transmission  Mueller  matrix 
corresponding  to  the  downward  propagation  by  Nd,  the  backscatter  Mueller  matrix  (also  known 
as  the  Stokes  reflection  matrix)  associated  with  the  point  target  by  M,  and  the  transmission 
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Mueller  matrix  corresponding  to  the  upward  propagation  by  Nu.  Suppose  the  radar  transmits 
a  wave  with  Stokes  vector  st.  After  passing  down  through  the  trees,  the  transmitted  Stokes 
vector  is  Nc&t,  and  after  bouncing  off  the  point  target,  the  scattered  Stokes  vector  is  (invoking 
Equation  3.12)  MNdst.  Finally,  after  passing  back  up  through  the  foliage,  the  scattered  Stokes 
vector  is  NuMNdat.  If  we  use  the  Stokes  vector  sr  to  describe  the  receive  antenna,  then  the 
received  power  is 


P 


where  D 


-sjDa, 

2 

NuMNd 


(4.2) 


The  objective  of  the  foliage  depolarization  measurement  process  is  to  estimate  Nu  and  Nd 
from  a  set  of  measurements  on  point  targets  located  under  the  foliage.  To  show  how  this  is 
done,  we  first  develop  some  mathematical  tools,  and  then  we  present  an  algorithm  for  inferring 
Nd  given  Stokes  vector  measurements  of  four  particular  PARCs. 


4.2  Mathematical  Tools 

In  this  section  we  will  develop  some  mathematical  tools  that  will  be  used  in  Section  4.3 
to  develop  the  measurement  algorithm.  We  will  derive  the  Stokes  reflection  matrix  (i.e.,  the 
backscatter  Mueller  matrix)  for  a  PARC,  and  then  will  show  that  the  Mueller  matrices  describing 
depolarization  during  downward  and  upward  foliage  penetration  are  related  by  Nd  =  N*  (so 
that  the  foliage  depolarization  measurement  problem  reduces  to  that  of  finding  only  one  or  the 
other  of  Nu  and  Nd). 


4.2.1  Stokes  Reflection  Matrix  for  a  PARC 

A  PARC  is  an  active  device  that  receives  a  radar  signal  with  one  antenna  (usually  linearly 
polarized),  amplifies  the  waveform,  and  retransmits  through  a  second  antenna.  The  polarization 
scattering  matrix  of  this  device  is  singular,  because  the  polarization  form  of  the  retransmitted 
field  is  independent  of  the  incident  field.  The  scattering  matrix  can  thus  be  represented  by  the 
outer  product  of  the  complex  two-vectors  representing  the  transmit  and  receive  polarizations: 


A 


t  r‘ 


(4.3) 
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where  the  tilde  indicates  a  vector  expressed  in  a  left-handed  coordinate  system. 

The  backscatter  Mueller  matrix  (Stokes  reflection  matrix)  of  the  PARC  can  be  found  using 
results  from  Section  3.3  plus  the  relations  T-1  =  (l/2)Tt  and  (a  <S>  b)T  =  aT  ®  bT: 


M  =  T'A  ®  A'T~X 

=  T‘{i  rT)®{tr})T-' 

=  ^ St{T'r®r')T 


No  matter  what  Stokes  vector  is  received  by  the  PARC,  the  transmitted  Stokes  vector  will 
be  of  the  form  i(.  To  see  this,  imagine  that  Stokes  vector  sT  is  incident  on  the  PARC.  The 
scattered  Stokes  vector  is  then  Msx  =  [(l/2)iti^]sx  =  [(1  /2)jjsx]it,  which  is  a  scalar  times 
the  vector  it. 

This  result  makes  sense:  the  PARC  always  transmits  the  polarization  st.  Note  that  the  Stokes 
reflection  matrix  is  singular,  as  is  the  polarization  scattering  matrix. 


4.2.2  Relationship  Between  Mueller  Matrices  for  Downward  and  Up¬ 
ward  Foliage  Penetration 

In  this  section  we  will  show  that  7Vd  =  Nj.  This  means  that  characterizing  foliage  depo¬ 
larization  requires  finding  only  one  or  the  other  of  these  Mueller  matrices. 

Consider  the  foliage  system  shown  in  Figure  4- la,  with  field  incident  from  above,  field  £* 
transmitted  down  through  the  foliage  to  the  area  beneath,  field  tc  incident  on  the  same  foliage 
from  below  after  having  bounced  off  a  target  under  the  foliage,  and  field  transmitted  up 
through  the  foliage.  The  latter  two  fields  are  expressed  in  a  left-handed  coordinate  system, 
consistent  with  the  monostatic  convention  described  in  Section  3. 

Figure  4- lb  shows  an  equivalent  four-port  network  representation  of  this  system,  with  Ports 
1  and  2  representing  the  horizontally  (H)  and  vertically  (V)  polarized  electric  fields  on  the  top 
side  of  the  foliage,  and  Ports  3  and  4  representing  the  H  and  V  fields  on  the  bottom  side  of 
the  foliage.  The  N-port  scattering  matrix  5  for  this  four-port  system  is  a  complex,  symmetric, 
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Figure  4-1.  Foliage  penetration  geometry,  and  equivalent  four-port  network. 
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4-by-4  matrix  relating  scattered  to  incident  waves  in  the  four  ports.  Let  us  express  S  in  terms 
of  its  two-by-two  submatrices,  denoted  here  by  the  dummy  matrices  P,  Q,  and  R: 


P  QT' 
Q  R 


(4.5) 


The  H  and  V  components  of  the  field  are  are  incident  waves  in  Ports  1  and  2,  while  the 
H  and  V  components  of  e*  are  outgoing  waves  in  Ports  3  and  4.  The  coupling  from  Pons  1 
and  2  to  Ports  3  and  4  is  described  by  the  lower  left  submatrix  Q  of  the  scattering  matrix  S : 


=  Qea  (4.6) 

Similarly,  the  H  and  V  components  of  the  field  ^  are  are  incident  waves  in  Ports  3  and  4, 
while  the  H  and  V  components  of  ed  are  outgoing  waves  in  Ports  1  and  2.  The  coupling  from 
Ports  3  and  4  to  Ports  1  and  2  is  described  by  the  upper  right  submatrix  QT  of  the  scattering 
matrix  5: 


U  =  QT£c  (4-7) 

Now  let  us  repeat  the  development  of  Section  3.3,  this  time  working  with  forward  scatter 
rather  than  backscatter.  We  will  extend  die  development  by  letting  the  foliage  transmission 
matrix  Q  be  random,  and  then  averaging  over  the  ensemble  of  all  the  random  matrices  Q  in  the 
probability  space  (this  ensemble  average  is  denoted  by  an  oveibar).  The  downward-transmitted 
Stokes  vector  is 


Ab  = 

=  TQ& Z  ®  Q'<& 

=  T{Q  ®  <?*)(£*  ®£;) 

=  TQ  ®  Q- 

=  ( \t'QWQ:T ')*. 

This  result  can  be  written  &  =  where 


(4.8) 
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(4.9) 


Nd  =  TQ^Q-T-1  =  ^TQ®Q'T' 

is  the  Mueller  matrix  relating  downward-transmitted  to  incident  Stokes  vectors.  Note  that 
both  the  incident  and  the  transmitted  Stokes  vectors  are  expressed  in  a  right-handed  coordinate 
system. 

Similarly,  the  incident  and  transmitted  Stokes  vectors  for  the  upward-traveling  waves  (Sc 
and  sd)  are  related  (as  expressed  in  a  left-handed  coordinate  system)  by 


id  =  T’ed  ®  £ 

=  T'QTec®Q^ 

=  T‘(QT®Q')(ic®e:) 

—  T~QT  ®  Qf  ec  ®  e'c 

=  {T‘QT 

=  (\T-QfTQ'TT)l  (4.10) 

This  result  can  be  written  £d  =  Nusc,  where 


Nu  =  T‘QT  ®  Ql(T')-1  =  ~T‘QT  ®  Q'TT  (4.11) 

is  the  Mueller  matrix  relating  transmitted  to  incident  Stokes  vectors  for  the  upward-traveling 
wave. 

The  matrices  Nu  and  Nd  can  now  be  related  to  one  another  by  direct  calculation: 


JVj  =  '-T(Qr«Q<fT< 

=  TQ®&(l-T<) 

=  T-Q  ®  Q'T~l 

=  Nd  (4.12) 
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We  now  have  Nd  in  terms  of  Nu.  Combining  Equations  4.2  and  4.12,  the  net  backscatter 
Mueller  matrix  for  the  foliage/target/foliage  combination  is 


D  =  NuMNd 
=  NjMNd 

=  NUMN*  (4.13) 

The  results  shown  in  Equations  4.4,  4.10,  and  4.13  will  be  used  in  Section  4.3. 


4.3  Measuring  Foliage  Depolarization  using  PARCs 

In  order  to  use  the  results  derived  in  Section  4.2,  let  us  suppose  that  a  radar  and  four  PARCs 
are  set  up  as  follows: 

1.  The  airborne  radar  transmits  a  single,  fixed  polarization  (say,  vertical  polar¬ 
ization).  Let  this  polarization  correspond  to  Stokes  vector  ^ .  This  Stokes 
vector  is  not  normalized,  i.e.,  it  has  a  first  element  which  is  not  equal  to 
unity,  but  rather  depends  on  the  (known)  radar  tranmsit  power. 

2.  The  radar  illuminates  four  PARCs  that  have  been  placed  under  foliage. 

3.  The  PARCs  are  built  so  that  their  receivers  are  all  matched  to  the  radar 
transmit  polarization  i0.  The  PARC  receiver  gain  will  be  described  by  the 
constant  of  proportionality  /?,  so  that  the  receive  Stokes  vector  common  to 
all  of  the  PARCs  is  /3s0. 

4.  The  PARCs  have  four  different  transmit  polarizations,  identified  by  Stokes 
vectors  x,,i  =  1,  ...,4. 

With  these  assumptions,  the  net  backscatter  Mueller  matrix  for  the  foliage/PARC/foliage  is 
(for  the  ith  PARC): 


Di 


(4.14) 


and  the  backscattered  Stokes  vector  from  the  ith  PARC,  which  we  denote  by  is 
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Li  =  Aa o 

=  - Nuii(l3s0)TNLs0 

=  (4.15) 

where  a  =  ^{0so)T s0  is  an  unknown  scalar  multiplier  that  depends  on  attenuation  through 
the  foliage  to  the  PARC  receive  antenna. 

The  estimation  of  Nu  will  be  achieved  in  two  steps:  first,  the  estimation  of  a  normalized 
Nu  (Section  4.3.1),  and  then  the  estimation  of  the  absolute  attenuation  corresponding  to  a 
(Section  4.3.2).  Since  the  PARC  receive  antennas  are  all  matched  to  the  radar  transmit  antenna, 
enough  power  should  be  available  to  the  PARCs  to  result  in  a  signal-to-noise  ratio  adequate 
for  parameter  estimation 


4.3.1  Measurement  of  Relative  Depolarization 

Equation  4.15  is  really  four  equations,  one  for  each  PARC.  Let  us  combine  them  into  one 
matrix  equation.  To  do  this,  arrange  the  vectors  /.  to  form  matrix  F,  and  also  arrange  the 
vectors  xt  to  form  matrix  X.  Next,  define  the  matrix  N  —  aNu.  Then  the  four  equations  can 
be  expressed  compactly  as 


F  =  NX  (4.16) 

The  matrices  F  and  X  are  both  known:  the  X  matrix  simply  expresses  our  specifications 
for  the  four  PARCs,  and  the  F  matrix  is  measured  directly.  Thus,  we  can  solve  for  the  desired 
matrix  N: 


N  =  FX -1  (4.17) 

The  unknown  matrix  Nu  is  then  proportional  to  TV.  In  section  4.3.2,  we  will  discuss  how 
to  determine  the  constant  of  proportionality  a. 

Strictly  speaking,  in  order  for  this  scheme  to  work,  all  we  need  is  any  set  of  four  linearly 
independent  s.x  that  span  the  space  of  real  four-vectors  and  that  therefore  form  an  invertible 
matrix  X.  However,  to  ensure  that  the  matrix  X  is  well-conditioned  (and  therefore  has  an 
inverse  which  is  robust  in  the  presence  of  noise),  we  should  pick  four  that  are  as  different 
from  each  other  as  possible.  A  good  set  would  be  H,  V,  45  degree  linear  polarization,  and 
left-circular  polarization,  with  Stokes  vectors  that  can  be  found  from  [3]  as: 
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(4.18) 


A  word  on  the  choice  of  left-circular  polarization  is  appropriate  here.  It  is  necessary  to  put 
some  energy  into  the  last  element  (£4)4  of  the  fourth  transmit  Stokes  vector  £4  in  Equation 
4.18,  in  order  to  stimulate  the  last  column  of  Nu  in  Equation  4.15.  PARCs  that  transmit  linear 
polarization  (for  example,  xlt  £2,  and  £3  in  Equation  4.18)  have  no  energy  in  the  last  element 
of  their  Stokes  vectors,  and  are  thus  unable  to  stimulate  the  last  column  of  Nu. 


While  it  may  be  difficult  to  construct  a  circularly  polarized  patch  antenna  such  as  that 
needed  to  transmit  Stokes  vector  £,,  it  is  not  impossible  [7],  and  at  least  one  PARC  vendor 
has  expressed  a  willingness  to  build  a  PARC  that  uses  this  type  of  antenna. 

This  choice  of  PARCs  gives  us  an  X  matrix 


I' 
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with  an  inverse 
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To  summarize,  the  depolarization  measurement  approach  proposed  here  involves  the  use  of 
four  PARCs,  all  of  which  are  designed  to  receive  vertical  polarization  but  which  transmit  H, 
V,  45  degree  linear,  and  left-circular,  respectively.  The  first  three  are  standard-issue  PARCs 
(VH,  VV,  V/45-degree),  while  the  last  (V/left-circular)  may  involve  some  effort  in  modifying 
or  redesigning  an  existing  PARC. 


4.3.2  Absolute  (Radiometric)  Measurement 

In  the  preceding  section,  we  showed  how  to  estimate  the  one-way  Mueller  matrix  of  foliage 
up  to  but  not  including  an  arbitrary  scalar  multiplier  a:  after  performing  this  partial  calibration, 
we  know  aNu  but  not  Nu.  To  determine  a  (and  therefore  Nu),  the  following  steps  would  be 
taken: 
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1 .  Measure  the  Mueller  matrix  Dx  for  the  W  PARC. 

2.  Remove  the  effect  of  relative  depolarization  by  pre-  and  post-multiplying 
Di  by  the  inverses  of  the  normalized  transmission  Mueller  matrices  deter¬ 
mined  as  shown  in  Section  4.3.1. 

3.  Complete  the  determination  of  a  by  using  the  known  magnitude  of  the  W 
PARC  response. 

The  Mueller  matrix  for  the  W  PARC  is 


A  =  K  (4.19) 

We  adjust  Dx  to  remove  depolarization  effects  by  premultiplying  Dx  by  the  inverse  of  N  = 
aNu  and  postmultiplying  by  the  inverse  of  NT.  The  result  is  the  relative-polarimetric-adjusted 
version  D 


D,  =  (aN„)-'D,(aNZ)-' 

=  i,(/51o)r  (4-20) 

But  we  know  the  values  of  the  vectors  on  the  right  side  (they  are  the  transmit  and  receive 
Stokes  vectors  for  the  W  PARC,  which  is  itself  calibrated),  and  we  can  measure  the  matrix 
on  the  left  side;  thus  we  can  solve  for  the  attenuation  a.  Let  the  first  elements  of  the  vectors 
on  the  right  hand  side  be  (ij)!  and  (/3iJ)i,  respectively,  and  the  upper  left  element  of  the 
matrix  Dx  be  (£>i)n-  Then 


(A)n  =  ^(ii)i(/3io)i 


a  = 


\ 


UiMgg) 

2(A) 


11 


(4.21) 


We  have  now  completed  the  measurement  of  the  upward-traveling  Mueller  matrix,  since  we 
know  aNu  and  a,  and  can  therefore  solve  for  Nu.  As  was  shown  in  Section  4.2.2,  this  also 
gives  us  Nj. 

The  measurement  of  depolarization  during  foliage  penetration  is  now  complete,  since  both 
the  upward-traveling  and  downward-traveling  Mueller  matrices  have  been  found. 
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APPENDIX  A 

STOKES  VECTORS  FOR  RIGHT-  AND  LEFT-HANDED 

SYSTEMS 


The  Stokes  vector  s  is  defined  in  [3]  in  terms  of  ex  and  ey,  the  elements  of  the  complex 
two-vector  describing  the  electric  field  (expressed  in  a  right-handed  coordinate  system)  by: 


So 

'  ie*|2  +  1  ew  1 2  ' 

Si 

\ex\2  -  |eu|2 

S2 

2Re[eIe’] 

s 3 

_  — 2Im[ese‘l  _ 

Section  8.2  of  [3]  (using  an  idea  from  [1])  shows  how  this  expression  can  be  more  elegantly 
expressed  using  a  Kronecker  product: 


5 

where  E_ 


and  where  T 


T(E®E-) 

fe-1 


10  0  1 

10  0-1 
0  110 
0  j  -j  0 


(A.2) 


Section  8.3  of  [3]  shows  that,  in  a  left-handed  system,  the  definition  of  s3  must  be  changed 


from  -Imfe,e* 


ization  must  be  1  0  0  1 


to  +Im[eKe*]  because,  for  example,  the  Stokes  vector  for  left-circular  polar- 

T 


regardless  of  whether  the  complex  two- vector  describing  the 
left-circular  field  is  expressed  in  a  right-handed  or  left-handed  coordinate  system  (the  repre¬ 
sentations  being  1  j  ]T  °r  1  —j  ]T,  respectively).  As  a  result  of  this  re-definition 
of  s3,  one  must  use  the  conjugate  of  T  instead  of  T  when  expressing  a  Stokes  vector  in  a 
left-handed  coordinate  system: 


I 

where  E 


T'{Z®R) 
.  ev  . 


and  where  T’  is  the  conjugate  of  the  matrix  T  defined  in  Equation  A.2. 


(A.3) 
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